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Abstract 

Using the graphical calculus and integration techniques introduced by the authors, we study 
the statistical properties of outputs of products of random quantum channels for entangled inputs. 
In particular, we revisit and generalize models of relevance for the recent counterexamples to the 
minimum output entropy additivity problems. Our main result is a classification of regimes for 
which the von Neumann entropy is lower on average than the elementary bounds that can be 
obtained with linear algebra techniques. 

1 Introduction 

As in classical computer science, randomized proofs and constructions are ubiquitous in quan- 
tum information. Since quantum mechanics is non commutative, the random objects of study 
are matrices. Therefore, quantum information theory provides a rich source of random matrix 
problems. 

One of the most important classes of problems in the mathematical aspects of quantum in- 
formation theory is the study of data transmission through noisy quantum channels. A famous 
conjecture reduced the calculation of the channel capacity for classical data to the question of the 
additivity of Minimum channel Output Entropy. The conjecture was stated in 1999 by King and 
Ruskai 1 and shown to be equivalent to the additivity of the Holevo capacity (and to other quan- 
tities of interest) by Shor [2J. For a long time, no counterexamples were available and additivity 
was proven to hold in many cases. A stronger, LP version of this question was also available and 
relevant to operator algebra and operator space theory. This version was disproved by Hayden 
and Winter in 2007, for all p > 1 [3J. The original conjecture, regarding von Neumann entropies, 
was disproved by Hastings in 2008 [4]. His very innovative argument exploited the idea of tubu- 
lar neighborhoods so as to considerably refine available estimates on random quantum channels. 
However, constructive, non-random counterexamples to any of these conjectures are still elusive. 

The random counterexamples to the various forms of the additivity conjecture rely on bounds 
on the Minimum Output Entropies (MOE) for single and product channels that follow mainly 
from two important ideas. Let <E> be a random quantum channel between matrix spaces such that 
the dimensions of the input and output spaces are large enough. The first key idea is that, with 
high probability, the Minimum Output Entropy of <i> is almost maximal: all output states are 
highly mixed. On the other hand, if one considers the product channel <& ® $ (where <E> is obtained 
by replacing the Stinespring unitary U defining the channel by its conjugate), then if one takes 
a maximally entangled (or Bell) state as an input, the output density matrix has always a large 
eigenvalue. This second important fact was observed by Winter, and it implies that the output 
state in question has low entropy, allowing for a violation of additivity. 
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Our work addresses bounds for conjugate product channels, and improves them in several cases. 
We provide a complete spectral description of the output of product channels when the input is 
maximally entangled. In [5] and [6], we have studied situations when the channels are conjugate 
((J> g) $) or independent ("1> g) VP) in two different asymptotic regimes. 

In this work, after recalling the aforementioned results and reviewing the techniques used in 
deriving them, we consider more general models of random quantum channels, from two different 
perspectives. We first generalize the linear scaling asymptotic regime to include the situations 
where the dimension of the input space is different from the dimension of the output of a quantum 
channel; however, all three parameters (the respective dimensions of the input, output and ancilla 
spaces) scale linearly. Then we move beyond the linear regime, considering situations where the 
dimensions of the output space and of the ancilla space scale in a non-linear fashion. Motivated 
by the search of improved bounds one may use in the study of additivity questions, we compute 
asymptotic expressions for the von Neumann entropies of output matrices for the models under 
consideration. 

The paper is organized as follows: in Section [5] we first review the tools available to study 
moments of outputs of random quantum channels. These techniques were introduced in [7J [8] and 
[5] and their first applications to quantum information theory were developed in [5J E] and [9]. In 
Section [3] we generalize the results of [5] to the case where the relative dimensions of the input 
and the output are different. In Section [4] we generalize the setting of [6] to the case where the 
relative dimensions of the input space and the ancilla space have relative polynomial growth. This 
is motivated by the recent results of [TO], where the authors consider the case k ~ n 1//2 (n being 
the dimension of the input/output space and k being the dimension of the ancilla space). We show 
that depending on the growth, different results occur and that the case where the ancilla space 
and the input space have the same dimension has a potential for yielding a bigger violation for 
the additivity of the entropy. Finally, in Section [5] we use these results to provide new bounds for 
von Neumann entropy of the output of product random quantum channels. 

2 Studying moments of outputs of random quantum 
channels: techniques and first examples 

In this section, we recall, for the convenience of the reader and for the sake of being self-contained, 
techniques to compute the eigenvalue distribution of random quantum channels, as well as a few 
results obtained recently with these techniques. 

The techniques rely on Weingarten calculus (subsection 12. 1[) and on a graphical model (sub- 
section [2]2} . Then, in subsections 12.41 and 12.51 we recall two applications of these techniques. 

2.1 Weingarten calculus 

In this section, we recall a few facts about the Weingarten calculus, useful to evaluate averages 
with respect to the Haar measure on the unitary group. 

Definition 2.1. The unitary Weingarten function Wg(n, a) is a function of a dimension param- 
eter n and of a permutation a in the symmetric group S p on p elements, defined as the pseudo- 
inverse of the function o v+ n* CT under the convolution for the symmetric group (if a denotes the 
number of cycles of the permutation a). 

Note that the function a *-> n* a is invertible for n ^ p, (to see that it is invertible for n large 
enough, observe that it behaves like n p S e as n —¥ oo). In this case, we can replace the pseudo- 
inverse by the inverse. We refer to [S] for historical references and further details. We shall use 
the shorthand notation Wg(u) = Wg(n, a) when the dimension parameter n is obvious. 

The following theorem relates integrals with respect to the Haar measure on the unitary group 
U(n) and the Weingarten function Wg. (see for example [11]): 

Theorem 2.2. Let n be a positive integer and (ii, . . . ,i p ), (i[, . . . ,i' p ), (ji, ■ ■ ■ ,jp), {ji, ■ ■ ■ ,j' p ) be 
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p-tuples of positive integers from {1, 2, . . . , n}. Then 



f u iljl --.u ipjp u i , ijr ..u i , j , dU = 

JU(n) 

8i.ii ...Si it 8i.ii ...8i ii Wgfn,™ 1 ). (1) 

£T,TG5 p 

If p ^ p' then 

( u iljl ---u ipjp u^7---W~r du = o. (2) 

J»(n) p p 

We are interested in the values of the Weingarten function in the limit n — > oo. The following 
result encloses all the data we need for our computations about the asymptotics of the Wg function; 
see [11] for a proof. 

Theorem 2.3. For a permutation a EE S v , let Cycles(cr) denote the set of cycles of a. Then 

Wg(n,a) = (-ir-# CT Wg(n, C )(l + 0(n- 2 )) (3) 

cGCyclcs(cr) 

and 

Wg(n,(l,...,d)) = (-l) d - 1 Cd _ 1 in-j)- 1 (4) 

-d+l^j^d-l 

where Ci — pjjrnj ,i * s the i-th Catalan number. 

The Catalan numbers and Wg are related to the Moebius function on the lattice of non-crossing 
partitions, as follows: 

Wg(n,a) = n- (p+H) (Mob(o-) + 0(n- 2 )) (5) 

where \a\ = p — #<r is the length of a, i.e. the minimal number of transpositions that multiply to 
a. We refer to [8] for details about the function Mob. 



2.2 Planar expansion 

The purpose of the graphical calculus introduced in [5] is to yield an effective method to evaluate 
the expectation of random tensors with respect to the Haar measure on a unitary group. In 
graphical language, a tensor corresponds to a box, and an appropriate Hilbertian structure yields 
a correspondence between boxes and tensors. However, the calculus yielding expectations only 
relies on diagrammatic operations. 

Each box B is represented as a rectangle with decorations on its boundary. The decorations 



are either white or black, and belong to S(B) U S*(B). Figure 1(a) depicts an example of boxes 
and diagrams. 

It is possible to construct new boxes out of old ones by formal algebraic operations such as sums 
or products. We call diagram a picture consisting in boxes and wires according to the following 
rule: a wire may link a white decoration in S(B) to its black counterpart in S*(B). A diagram 
can be turned into a box by choosing an orientation and a starting point. 

Regarding the Hilbertian structure, wires correspond to tensor contractions. There exists an 
involution for boxes and diagrams. It is antilinear and it turns a decoration in S(B) into its 
counterpart in S*(B). Our conventions are close to those of |121 113]. They should be familiar 
to the reader acquainted with existing graphical calculus of various types (planar algebra theory, 
Feynman diagrams theory, traced category theory). Our notations are designed to fit well to the 
problem of computing expectations, as shown in the next section. In Figure [l(b)| |l(cj| and [T(d)] 
we depict the trace of a matrix, multiplication of tensors and the tensor product operation. For 
details, we refer to [5]. 

The main application of our calculus is to compute expectation of diagrams where some boxes 
represent random matrices (e.g. Haar distributed or Gaussian). For this, we need a concept of 
removal of boxes U and U. A removal r is a way to pair decorations of the U and U boxes 
appearing in a diagram. It therefore consists in a pairing a of the white decorations of U boxes 
with the white decorations of U boxes, together with a pairing j3 between the black decorations of 
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Figure 1: Basic diagrams and axioms: (a) diagram for a general tensor M; (b) trace of 
a (l,l)-tensor (matrix) M; (c) Scalar product (y | M \ x); (d) tensor product of two 
diagrams. The labels round, square and diamond-shaped labels correspond to pairs of 
dual finite dimensional complex Hilbert spaces. 



U boxes and the black decorations of U boxes. Assuming that T> contains p boxes of type U and 
that the boxes U (resp. U) are labeled from 1 to p, then r = (a, P) where a, j3 are permutations 
of S p . 

Given a removal r 6 Rem(X>), we construct a new diagram T> r associated to r, which has the 
important property that it no longer contains boxes of type U or U. One starts by erasing the 
boxes U and U but keeps the decorations attached to them. Assuming that one has labeled the 
erased boxes U and U with integers from {1, . . . ,p}, one connects all the (inner parts of the) white 
decorations of the i-th erased U box with the corresponding (inner parts of the) white decorations 
of the a(i)-th erased U box. In a similar manner, one uses the permutation ft to connect black 
decorations. 

In [5], we proved the following result: 

Theorem 2.4. The following holds true: 

Ev{T>)= ©rWgKa/T 1 ). 

r=(a,3)eR,emtHZ>) 

2.3 Wishart matrices, Marchenko-Pastur distributions and their 
entropy 

We recall the definition of a free Poisson (or Marchenko-Pastur) random variable [14] . For c > 0, 
the probability measure 

i/4c — (x — 1 — c) 2 
tt c = max(l - c,0)<5o + — l[i+c-2^,i+ c +2^](z) dx 

is called a free Poisson measure of parameter c. The plots of the densities for these measures are 
plotted in Figure [2] 

The free Poisson distribution arises in random matrix theory as the almost sure limit of the 
eigenvalue counting distribution for Wishart matrices, i.e. matrices X n X* where X n is an n x [cn\ 
matrix whose entries are i.i.d. standard complex Gaussian random variables of variance 1/n. 

From a combinatorial perspective, the free Poisson distribution n c has the nice property that 
all its free cumulants are equal to c. Hence, the free moment-cumulant formula (see [15], Lecture 
11, pp. 173) reads 

I x p d-K c (x) = ]T c #CT , (6) 

where #<r denotes the number of blocks of the non-crossing partition a. From the moment formula, 
one can obtain the value of the following integral, useful in the computation of von Neumann 
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Figure 2: Densities for the Marchenko-Pastur measures of parameters c— 1/5, c — 1 and 
c = 5. For c = 1/5, only the absolutely continuous part of the measure was plotted; 7Ti/5 
has a Dirac mass of 4/5 at x — which is not represented. 



entropies [IS] : 
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ifO<c< 1. 
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2.4 Application 1: Fixed ancilla space 

The counterexamples to the additivity conjecture obtained so far arise from the random choice of 
a quantum channel from the ensemble 



$ : Mr 



Mr. 



(8) 



given by $(X) = Tr k [U(X ® |0}(0|)t/*], where U € bt(nk) is a random unitary matrix and |0}(0| 
is an ancilla rank-one projector. The counterexamples rely on the idea of using a tensor product 
of conjugate channels <E> ® $ and, more precisely, the output of this channel when the input is the 
Bell state ({ei, . . . , e m } is some fixed basis of C m ): 



E r , 



1 m 



|e«><ej 



(9) 



In the regime where m = infc with n and ft integers and £ £ (0, 1) a fixed parameter, we gave in 
[5] a complete spectral description of the (random) density matrix [$ ® "3>](B tn fc). The following 
result improves on the previously known bound of Winter, Ai ^ t: 

Theorem 2.5. Almost surely, as n — > oo, the k 2 non-zero eigenvalues of the random matrix 
[3> ® $\(E t nk) £ M n 2(C) converge towards the deterministic probability vector 
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A^n, corresponding to kt — 1, the eigenvalues are t + -p- 
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In the particular case $ : A'f 
with multiplicity one and tj — tj, with multiplicity k 2 " — 1. Not only the value of the largest 
eigenvalue is improved from ^ to x + — t^! but the lower eigenvalues are also computed. 

A better understanding of the spectrum of the output matrix for the product channel yields 
immediately better bounds for the Minimum Output Entropy of $ ® <E>. Applications of this result 
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are twofold. First, it allows for violations of additivity of Renyi entropies, for all p > 1, just by 
using a qubit as an ancilla state for the input (t = 1/2). Second, the improvements on the bound 
for the entropy of the product channel can yield better minimum values of k needed to obtain 
violations of additivity for p — 1 (see [171118] ). 



2.5 Application 2: Ancilla and input space of linear dimensions 



Theorem 12.51 is highly non- intuitive: it is not clear why the small eigenvalues should all behave in 
the same way. Moreover, there was numerical evidence [5] that the spectrum should not be flat 
beyond the second eigenvalue. This raises the question of what happens when the ancilla space is 
not of fixed dimension, but rather of dimension comparable to the input space. 

The study of such asymptotic regimes was initiated in [BJ. After stating the main result obtained 
in that paper, we shall generalize in Section [3] the setting, allowing for the dimension of the input 
space to vary linearly with the dimension of the output. 

In [BJ, Section 6.3, we considered random quantum channels $ : M n (C) —¥ M n (C) obtained 
from random Haar unitary matrices U G U(nk) via the Stinespring representation 



$(X) =Tr fc \U{X®P k )U* 



(10) 



where Pk G Mk(C) is a non-random rank-one projector (pure state) and both n and k grow to 
infinity at a constant ration k/n — > c > 0. The diagram for such a channel is represented in Figure 
® 



c $(X) • 




Figure 3: Diagram for a quantum channel with equal input and output spaces. The 
state Pk of the ancilla space is omitted, since it has no role to play in the computations. 
Round labels attached to boxes correspond to input/output spaces C™ and square symbols 
correspond to ancilla spaces C k . 



For the regime where both n and k grow to infinity at a constant ratio c, the main result of [BJ 
is as follows. 

Theorem 2.6. Consider a pair of conjugate random quantum channels <E>, <E> in the regime where 
n,k —¥ co, k ~ cn. The eigenvalues \\ ■ ■ • ^ A n 2 of the random matrix Z n — [$ ® <&](I? n ) are 
such that: 

• The first eigenvalue satisfies cn\\ — » 1 (in probability). 

• The distribution 53™_ 2 8 c 2 n 2 X . converges a.s. to a free Poisson distribution of parameter 
c 2 . 

Here we see that a new phenomenon of two different convergence rates for eigenvalues appears. 
This is due to the fact that the U — U model for the product channel contains the "conjugation" 
symmetry. Instead of considering a channel and its complex conjugate, we looked in [6j at two 
independent quantum channels, taken from the same ensemble. 

Theorem 2.7. In the regime k ~ cn, n —¥ co, let Z n — [$ ® ^ , ](£' ?l ) be the output of the product 
of two independent quantum channels $ and ty, when the input is a maximally entangled state E n . 
Then, almost surely, the distribution of the rescaled output matrix c n Z n converges towards a free 
Poisson of parameter c 2 . 

A striking feature of this asymptotic regime is that the von Neumann entropies of the U — V 
and U — U models are almost the same. It is then natural to ask whether the U — U symmetry is in 
fact needed to obtain violations of the additivity. Indeed, it seems that the largest eigenvalue for 
the output of the product channel does not play a big role in the bounds for the entropy. Having 
counterexamples with independent channels will be an important step to a better understanding of 
additivity violations. Possible violations with independent channels would be generic, as opposed 
to conjugate-channels violations which are not. Also, these considerations give concrete hope that 
larger violations of additivity could be achieved. 
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3 Generalized linear setting — input and output spaces 
of different dimension 



In this section, we generalize the model of quantum channels we have considered, by allowing input 
and output spaces of different dimensions. We consider random quantum channels <£> : A4 m (C) — > 
M n (C) defined by 

$(X)=Tr k [U(X®Pi)U*], (11) 
where I — nk/m and Pi is a deterministic rank-one projector in Mi(C); we tacitly assume that 
I is an integer. All three dimensions m,n and k grow to infinity, at constant ratios: m/n — s> b 
and k/n — > c, with 6, c £ (0, oo) fixed constants. The generalized diagram corresponding to $ ias 
depicted in Figured] 
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Figure 4: Diagram for a quantum channel with different input and output spaces. Round 
labels attached to boxes correspond to output spaces C™, square symbols correspond to 
ancilla spaces C fc and diamonds correspond to input spaces C m . The rank-one projector 
Pi is omitted. 



When presented with the maximally entangled (or Bell state) E m G Mm 2 (C) as an input, the 
product conjugate channel $ ® $ produces a random density matrix 



Z= [$®$](£ m ) eM n 2(C). 



(12) 



The remaining of this section is dedicated to the study of the random matrix Z, depicted in 
Figure [5] The analysis of the spectrum of Z follows closely corresponding results in [6] , Section 
6.3, which is a specialization of this section, in the case m — n (i.e. 6=1). The spectral properties 
of the output random matrix Z are summarized in Theorem 13.41 the main result of this section. 
The reader in invited to compare the conclusions of Theorems l2.6l and l3.4l 
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Figure 5: Diagram for the output of a product of two conjugate channels, when the input 
is the maximally entangled state. The complex Hilbert spaces associated to labels are as 
follows: o • C", □ ■ ~» C fe , o ♦ ~» C m and a a ~> C l . 



The first step of our analysis of the output matrix is the computation of the asymptotic 
moments. 

Proposition 3.1. Consider a sequence of random quantum channels $ n where m,n,k — > oo, 
m/n — > b and k/n — >■ c. The asymptotic moments of the output matrix Z — [$ (g> $](_E m ) are given 
by: 



ETr 



Tr(Z) = 1 



1 c 



■ Tr l('^) P ] = l + o(l), Vp^ 3. 
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Proof. The starting point of the proof is the following exact formula for the moments of the random 
matrix Z, obtained from the graph expansion technique detailed in [5]: 

E[Tt{Z p )} = fc #a n #(a7_1) m #(,3i) ^Wg(a/?" 1 ), (13) 

a,/365 2p 

where the permutations y, 5 € & P are defined as follows. We relabel the index set {1, 2, . . . 2p} 
as {p B , ■ ■ ■ , I s , 1 T , 2 T , . . . ,p T } in order to make precise the association of indices with blocks 
corresponding to the "top" channel <E> (l T ,...,p T ) and blocks corresponding to the "bottom" 
channel <!> (I s , . . . ,p B ). With this notation, we define the permutations 

7 (i T ) = {i-l) T , y{i B ) = (i + l) B and 5{i T ) = i B , 5(i B ) = i T . (14) 

Using the asymptotic expressions for the dimensions m ~ bn, k ~ an and for the Weingarten 
function (see Theorem 12. 3|) 

Wgia/3- 1 ) ~ {nky 2p -^^ Mob(a/3~ 1 ) ~ c -^P-\^- x \ n -^-n^\ Uoh ( a p-^ ( 15 ) 

equation (|13[) becomes 

E[Tr(Z p )]~ ^2 c~ (|a| + |Q ' 9_1|) 6 p - |,35| n- /(a ' /3) Mob(Q^" 1 ), (16) 

where the function f(a, p) is given by 

f(a,0) = \a\ + \aj- 1 \ + \pS\+2\ a p- 1 \ -p. (17) 

In order to find the dominating terms in the sums (|13p or (|16|1 , one has to minimize the quantity 
f(a, P) over the permutation group S2 P - This has been done in the proof of Theorem 6.8 of [6]: 

• for p — 1, f(a, p) ^ 0, with equality iff. a — P = id; 

• for p = 2, /(a, P) ^ p, with equality iff. a — P £ {id, 5,7}; 

• for p ^ 3, f(a, P) ^ p, with equality iff. a — P — 8. 

One concludes now by plugging the optimal values for the permutations a and P into equation 
Gil). □ 

Theorem 13 . 1 1 only gives a partial description of the spectrum of the random matrix Z; from the 
moment information one can deduce that there are some eigenvalues on the scale n _1 and that the 
rest of the spectrum is distributed on lower scales, such as n~ 2 . Obtaining information on the lower 
scale eigenvalues by brute force via the moment method is a difficult task, since their asymptotic 
contribution is negligible with respect to the contribution of the eigenvalue(s) on the scale n _1 . 
The trick we using to obtain information about the smaller eigenvalues is inspired by Hayden and 
Winter's proof of the existence of a large eigenvalue. Their proof contains, as a byproduct, some 
information on the eigenvector for the large eigenvalue. We introduce the orthogonal projection 
Q = I —E, where E = E n is the maximally entangled state on a product of two copies of the 
output space C™ <g> C™. Using the (rank n 2 — 1) projector Q, we shall obtain some information on 
the smallest n 2 — 1 eigenvalues of the random output matrix Z, by analyzing the "compressed" 
matrix QZQ (which, in a suitable basis, corresponds to considering the (1, 1) minor of Q). 

Proposition 3.2. Almost surely, the matrix c 2 n 2 QZQ converges in distribution, to a free Poisson 
( or Marchenko-Pastur) law of parameter c 2 . 

Remark 3.3. By 'almost surely', we mean with probability one, in any probability space on which 
the whole sequence (indexed by the input dimension) of random quantum channels is defined. In 
this paper, we supply no proofs of almost sure convergence results, as they require further -not so 
enlightening- technicalities. We refer the interested reader to the appendix of JjS|/ for details. Let 
us just mention that the proofs rely on the Borel-Cantelli lemma. More precisely, one proves that 
the covariance of any moment behaves as 0{n~ 2 ) as the dimension goes to infinity. The fact that 
0(n~ 2 ) is summable over n makes it possible to use the Borel-Cantelli lemma. 
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To prove this result, one can use the arguments of Theorem 6.9 from [6j mutatis mutandis. 
The method of moments is employed by computing the moments of the random matrix c 2 n 2 QZQ 
and showing that they converge to the corresponding moments of the free Poisson distribution of 
parameter <?: 

lim \E[Ti{cnQZQY] = f x p dn c 2(x). (18) 

n— »oo Ti J 

The reader can note that the parameter b (describing the size of the entangled input) has 
no influence on the lower part of the spectrum. It only appears in the expression of the largest 
eigenvalue of Z, as stated in the following Theorem, which is the main result of this section. 

Theorem 3.4. Consider a pair of conjugate random quantum channels $ in the regime where 
m, n,k —¥ co, m ~ bn and k ~ cn. The eigenvalues Ai ^ • • • ^ A n 2 of the random matrix 
Z n — [$ ® $](Em) are such that: 

• The first eigenvalue satisfies (c/b)n\\ —¥ 1 in probability. 

• The distribution -r^i X]"=2 ^c 2 n 2 \i converges almost surely to a free Poisson distribution of 
parameter c 2 . 

The proof of this result combines Theorems 13.11 13.21 and Cauchy's interlacing theorem f |19|. 
Corollary III. 1.5). 

Remark 3.5. The almost sure convergence argument described in Remark \3.3\ does not extend to 
the first item of Theorem\3.4\ as the covariances tend to zero but are not summable. 



4 Non-linear output dimension 

In this section we generalize the "linear" model of [B] in a different direction than we did in Section 
[3] We shall consider random quantum channels $ : M n (C) — > M n (C) defined by the Stinespring 
representation (|10p. where the dimension of the ancilla space C fe scales with n in a non- linear 
fashion: 

k ~ cn d , when n — s- oo. (19) 
Here, c > and d > are two real parameters of the model. As before, we are interested in the 
spectral properties of the random matrix Z = [$ ® $](£„), where E n is the rank-one projection 
on the maximally entangled state in the input space of the product channel, C™ ® C™. 

Before performing a detailed analysis of the spectrum of Z, let us make some observations on 
the role of the parameters c and d, as well as on several particular cases already treated in the 
literature. Generally, in such models, the parameter c will play the role of a scaling parameter in 
the limiting spectral distribution. This phenomenon can be observed in [20], Theorem 5 or in |BJ, 
Theorem 6.11. On the other hand, the exponent parameter d will have a more qualitative role to 
play, as it will decide the type of behavior of the spectrum of the random matrix Z. 

Several particular cases of this very general model of random quantum channels have already 
been studied in the literature. In [5], we studied the case where the dimension k of the ancillary 
system is fixed, which corresponds in our setting to the case d = (hence k = e). It has been shown 
(Section 12.41 Theorem 12. 5p that in this situation, the random matrix Z has one large eigenvalue 
(Ai = + fc~ 2 — fc -3 )) and that the rest of the spectrum is "flat": A2 = •• ■ = A fc 2 = k~ 2 — k~ 3 . 
In [6j, the case d = 1 was investigated; this corresponds to a coupling with an ancilla space of 
dimension k which scales as k ~ cn. The situation was rather different in this case: one large 
eigenvalue of size (cn) -1 was observed, and the lower spectrum was not flat, having a free Poisson 
7r c 2 shape. As a final remark, note that the model under study here is different than the one in 
Section [21 where inputs of arbitrary size were considered. 

The main result of this section is the following theorem, which classifies the spectral behavior 
of the output random matrix Z in terms of the parameter d. 

Proposition 4.1. The asymptotic moments of the random output matrix Z = [$ ® <E>](_E n ) are 
given by: 

1. Ifd = (see J^): 

K[Tr(Z P )] =(\ + ^-^y + (c 2 -l)(^-^Y + °W V P >2. (20) 
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2. Ifde (0,1): 



Ifd=l (see 



I Ifde (1,2): 



E[Tr(Z p )] ~ 2c- 2 n- 2d - 2fc~ 2 i/p = 2; 

E[Tr(Z p )] ~ c" p n" dp ~ £T P Vp ^ 3. 

E[Tr(Z p )] ~ (1 + 2c" 2 )n" 2 if p = 2; 

E[Tr(Z p )] ~ C - p n- dp ~ fc" p Vp ^ 3. 



C[Tr(Z p )] ~ n- (2p - 2) l /p< 2 



(21) 



(22) 



2-d' 



E[Tr(Z p )] ~ (l + c- p )n~ dp ~ (l + c- p )n- (2p - 2) ^ p =_l_; (23) 



i[Tr(Z p )] ~ C - p n- dp ~ fc" p i/p > 2 



2-d 
5. I/d > 2: 

E[Tr(Z p )] ~ n" (2p " 2) Vp > 2. (24) 

Proof. The starting point of the proof is the following exact moment formula, obtained via the 
graphical calculus introduced in [5] and the Weingarten formula: 

E[Tr(Z p )] = ^2 k*° , n* la ' y ~ 1)+ * m - p Wg{ap- 1 )- (25) 

After using the asymptotic expression for the ancillary dimension k ~ cn d and the Weingarten 
function 

Wg(a/3" 1 ) ~ ( nk y*P-\<*f>-^v-\<*^\ Mob(a/3~ 1 ) 

_ c -2p-| Q /3-M n -(d+i)(2p+|^- 1 |) Mob(a/3- 1 ), 

we obtain the following expression (which holds if the right-hand-side is non-zero): 



(26) 



E[Tr(Z p )] ~ J2 c- (|Q| + |a/3 1|) n- s(a ' /3) Mob(a/3~ 1 ), (27) 

where S^a, /3) is the exponent of n in the preceding sum: 

S(a,/3) = d\a\ + \a 1 - 1 \ + \06\ + (d + 1) la/3" 1 1 - p. (28) 

In this proof, since the permutation 5 is a product of disjoint transpositions, thus an involution, we 
shall use the fact that \a5\ — ct5~ 1 1 = a _1 <5|. In order to find the dominating terms in the sum 
(|27|l . one has to minimize the quantity S(a, /3) over <S 2 p . The solution to this problem is obtained 
in three steps. First, using the following triangular inequalities: 

M + ia/r 1 ! > 

we obtain that 

S(a,p) > d\p\ + |/3 7 _1 | + W\ -p=: Sitf), (30) 
where the inequality can be saturated if , e.g. a = f3. The minimization problem for the Si 
function, although much simpler that the two- variable problem for S, can be further simplified by 
replacing the two-cycle permutation j with the full-cycle: 

minimize S 2 (/3) = d|/3| + Ifi^ 1 1 + \/35\ - p, (31) 
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where 7 = (p T • ■ ■ 2 T 1 T I s 2 s ■ • -p s ). Since the permutation S is an element of the geodesic 
id — ¥ 7, the solution for the S2 minimization problem is easy to find using one or more of the 
following inequalities: 

101 >0; 

l/3| + l/?7 _1 | > |Tl = 2p - 1; (32) 
I/35I + I/37" 1 ! > \5T 1 \=P-1- 
The solution for the 5*2 problem, in terms of the value of the parameter d, is summarized in Table 

m 



d 


Minorant for S 2 {/3) 


Equality cases 





-1 


{*->7> 


(0,2) 


dp-1 


S 


2 


2p- 1 


{id 5} 


(2,oo) 


2p- 1 


id 



Table 1: Solution to the minimization problem of equation f|3 1 [) . 



Next, we move towards finding the minimum of Si(/3), defined in equation (|30[1 . The permu- 
tations 7 and 7 are at distance one: 

7 = 7'(1 S P T ), (33) 
hence the same holds for /37 _1 and /?7 _1 . We have thus Si(fi) = Sziff) ± 1 and, even more 
precisely, 

J S'2(/3) — 1 if I s and p T are in the same block of /37 _1 , 
\S 2 (I3) + 1 otherwise. 

In the same manner that is was argued in [6], p T and I s belong to the same block of 07 -1 
if and only if f3 sj 7 (the permutations being compared with the partial order relation on the 
corresponding non-crossing partitions). Analyzing the different equality cases in Table [T] and 
using the fact that the unique element of the geodesic set {id — » 5} which is smaller than 7 is 
P — id, we conclude that, for d ^ 2, Si(/3) ^ 2p — 2, with equality iff. /3 = id. The case d G [0, 2) is 
more intricate, since one cannot have Si(f3) — Sa(/3) — 1 and saturate at the same time the lower 
bound for S2. 

For d = 0, it was shown in [5] that Si(/3) ^ 0, with equality iff /3 G {(5, 7}. Choose d G (0, 1] 
and consider f3 G 52 P such that Si(f3) = S2 (/3) — 1. Then, since j3 ^ 7, one has |/?<5| ^ p and thus 

Si(p) = 8i{p) - 1 = d\(3\ + I/37- 1 ! + -p - 1 

= d(|/?| + I^T'l) + (1 - d)l/97 _1 | + \0S\ -p-1 (35) 
^ d{2p- 1) + (1 - d) - 1 = d(2p - 2). 

We conclude that, for d G (0, 1], p ^ 3 and /3 such that Si(fi) = S , 2 (/3) - 1, Si(/3) > Si (6) = dp. 
Thus, for d G (0, 1] and p 3, the minimum dp is attained only at the point ft = 8. For p = 2, an 
exhaustive search in 54 reveals that, for d G (0, 1), <Si(/3) ^ dp, with equality iff (3 = 5,7, and, for 
d = 1, Si(/3) > p, with equality iff /? = id, 5, 7. 

In the case d G (1, 2), the situation is different. We shall consider three cases: 

1. P £ {id — ¥ 7}. Since f3 is not a geodesic permutation, we have that |/3| + |/?7 _1 | ^ 2p + 1, 
|/3| + |/35| > p + 2, |/8| > 1 and \/3S\ ^ 1. It follows that 

Si(/3) ^ S 2 (/3) - 1 > 101 + |/?7 _1 | + (d- 1)(I0I + I0<5|) + (2 - d)|/35| - p-1 

(36) 

^ 2p + 1 + (d - l)(p + 2) + (2 - d) - p - 1 = dp + d > dp = Si (5). 
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2. ft 6 {id — > 7}. In this situation, Si(/3) = S , 2(/?) — 1 and thus (we use the fact that \/38\ p 
in this case) 

Si(ft) =S a (fi)-l> \P\ + IPT'l + (d- 1)|/9| + - p - 1 
> 2p - 1 + +p - p - 1 2p - 2 = Si (id), 

with equality iff. ft = id. 

3. (3 £ {id ->• 7}, j9 £ {id ->■ 7}. In this situation, Si(ft) = S 2 (ft) + 1 and thus 

Si(ft) = S 2 (/?) + 1 > |/?| + |0 7 _1 | + (d- 1)(|/3| + |/3«|) + (2 - d)|/35| -p+ 1 
> 2p - 1 + (d - l)p + (d - 2)0 - p + 1 > dp = Si (5), 

with equality iff. /3 — 8. 

Analyzing the three cases, we conclude that if p ^ 2/(2 — d), we have Si(ft) ^ 2p — 2, otherwise 
Si (ft) ^ rfPi the unique minimizers being respectively ft = id, ft = S and ft € {id, 8} for the 
interface case p = 2/(2 — d). The answer to the Si minimization problem is summarized in Table 



d 


Minorant for S\(fi) 


Equality cases 








{^7} 


(0,1) (p = 2) 


2d 


5,7 


1 (p = 2) 


2 


id, <5, 7 


(0,1] (p>3) 


dp 




(1,2) (p<j^- d ) 


2p-2 


id 


(i-2) (p=^g) 


2p-2 = dp 


id, 5 


(1,2) (p >5 ^) 


dp 


8 


[2,oo) 


2p-2 


id 



Table 2: Solution to the S\ minimization problem of equation (|30[) . 



In order to solve the initial problem for the two- variable function S(a,ft), one needs to notice 
that, for a fixed value of ft, the only possibility to saturate the inequalities (|29[) at the same time 
is a = ft. Since, for d > 0, both these inequalities are used to go from S to Si, we conclude that 
for strictly positive d, the solution for the S minimization problem can be found in Tabled with 
a — ft. The solution for d = has been entirely described in [5]. One concludes by replacing the 
values for the minimizing permutations into equation (|27|l . □ 

In the cases d = and d 2, the behavior of the eigenvalues can be easily deduced from 
the moment information, since one can identify in the formulas the moments of some probability 
distribution. In the case d = 0, as it was argued in [5], the larges eigenvalue of Z converges to 
- -+- 4j — % and the k 2 — 1 others converge to 4j — Jj. Such a behavior is typical for the model 
we study, with a spectrum containing one large eigenvalue and k' 2 — 1 (or n 2 — 1) identical smaller 
eigenvalues. The case d ^ 2 is somehow atypical: all the n 2 eigenvalues behave like n~ 2 . This 
is due to the fact that the "large" eigenvalue, which usually behaves as has no contribution 
asymptotically. We summarize these facts in the following proposition. 

Proposition 4.2. In the regime d — (which corresponds to considering a fixed ancilla dimension 
k = c), the eigenvalues of the random matrix Z are such that, almost surely, in the limit n — > 00, 

Ai -> - + \ \ and A 2 , . . . , A fc2 -> \ \. (39) 

c c 2 c s c 1 c A 
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In the regime d 2, the (rescaled) empirical spectral distribution of Z converges to the Dirac mass 
at 1, Si: 

1 " 2 

lim -j V5„2 Xi =8i. (40) 

As in Section O in order to understand fully the eigenvalues of the random matrix Z, we need 
to understand the lower part of the spectrum, in the remaining cases d £ (0, 1) and d £ (1, 2) 
(the case d — 1 being treated in [6]). This is done by "pinching" the matrix Z by the projector 
Q — I — E n , orthogonal to the maximally entangled state -En- 
Proposition 4.3. The matrix k 2 QZQ converges, in moments, to the Dirac mass at 1, Si. 

Proof. We follow the idea of the proof of Theorem 6.10 in [jjj, compute the moments of the rank 
k 2 matrix k 2 QZQ, and show that they converge to the corresponding moments of the limit law: 

lim ^-E[Tr(k 2 QZQ) p ] = 1. (41) 

n— too K 

After replacing Q = I — E and developing the product, we get 

-!-E[Tr(fc 2 QZQf] = fc 2p ~ 2 E[Tr(I -E)Z(I -E)Z ■ ■ ■ (I -E)Z] 
k 

= c 2p n 2p - 2 Y, (-l) l/_1(B)l n- |/_1(J5)l E[Tr/(l)Z/(2)Z.../(p)Z], (42) 
f£F p 

where T is a set of the 2 P choice functions / : {1,2, . . . , p} — > {l,E}. The factor 

n -\r\E)\ is 

due to the normalization of the Bell states E = E n . The moment E[Tr f(l)Zf(2)Z ■ ■ ■ f(p)Z] is 
computed via the graphical Weingarten calculus: 

E[Tr f(l)Zf(2)Z ■■■ f{p)Z] = ^ k# a n* (af ~~ 1 )+#<W-p Wg(a^ 1 ), (43) 

where / £ <5>2 P is the permutation associated to the choice function / £ T v describing the way / 
connects the different instances of the channel (the arithmetic operations of indices i should be 
understood modulo p): 

j\j(i-l) T if/« = I 
if f(i) = E, 

. B f\(i + 1) B if/(i + l) = I 
\i T if/(» + l) = B. 

Exactly as in the proof of Theorem 6.10 from [5], one has to understand the possible cancel- 
lations of high powers in n. In order to do this, we rewrite the non- asymptotic equation (|42p 
as 



^E[Ti{k 2 QZQ) p ] = Yl fc 4p - |Q| - 2 n 3p - |/3i| Wgf^" 1 ) ^ (-l) lr "- ! „ 



a,/3es 2p /e^p 

(44) 

As in [5], we can show that for all permutations a £ V, the sum over all choices / £ T v is exactly 
zero, where 

V = {a £ 52 P 3i £ {1, . . . ,p} s.t. cr(i T ) = i B or a(i s ) = i T } 

(45) 

= {cr £ iS2 P | f<5 has at least one fixed point}. 



Hence, 



±E[Tr(k 2 QZQ) P ] ~ £ ( _ 1) |/- 1 ( B )l c 2 ! ,-2-(|a| + |^- 1 |) n - S (^,/) Mob(a ^-i )) (46) 

/ej^p, a,/36S 2 p, a£V 
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where 

S(a,0,f) = \/3S\ +d\a\ + (d+l)\ar 1 \ + + \af~ 'l - p(2d + 1) + 2d. (47) 

Since a £ V, aS has no fixed point, and hence \a5\ ^ p. Using the facts that |a/3 _1 | + \/3S\ ^ \aS\, 
|a/3 _1 | > and \a\ + laf' 1 ] > |/|, we obtain that 

S(a,p,f) > \f- 1 (E)\ + d\f\ + (l-d)\af- 1 \-2dp + 2d, 

with equality if and only if a = /? € {id —>-/}. 

The number of cycles of / is easily shown to be: 



2 if / = I, 

otherwise. 



#f = { ,,-!,„„ i_ ( 4 8) 



hence, + > 2dp - 2d, with equality iff / = I. We can conclude that S(a,/3,f) ^ 0, 

with equality if and only ifa = /3 = I = 7- Replacing these values in equation (|46|l . we obtain the 
announced result. □ 

Using Cauchy's interlacing theorem ([H], Corollary III. 1.5) for the eigenvalues of QZQ and 
those of Z, we obtain the complete description of the spectrum of the random output matrix. 

Theorem 4.4. In the regime k ~ cn d , with d € (0, 1), the eigenvalues Ai ^ • ■ ■ ^ A n 2 o/ 
Z = [$ g> $](£„) satis/?/.- 

• /n probability, k\\ —¥ 1. 

• Almost surely, k -i_ 1 X^=2 ^fc 2 A; converges in distribution to the Dirac mass at 1, Si. 

• TTie remaining n 2 — fc 2 eigenvalues are null: A fe 2 +1 = • • • = A n 2 = 0. 

We finally treat the regime d £ (1,2), stating the results and underlining the (small) differences 
between the proofs in this case and the proofs of Proposition 14.31 and, respectively, of Theorem 

H3 



Proposition 4.5. In the regime d G (1,2), the matrix n 2 QZQ converges, in moments, to the 
Dirac mass at 1, Si. 



Proof. As in the proof of Proposition 14.31 we obtain 



E[Tr(n 2 QZQ) p ] ~ £ (-l)''"^'^'-'^ 'V 5 '"'^ Mob^" 1 ), (49) 

feT p , a,l3eS 2p , a?V 



where 

S(a, /3, f) = \PS\ +d\a\ + (d+ TJla/T 1 1 + \r\E)\ + {af^l - 3p + 2. (50) 
Using the same arguments as before, we obtain that 

S(a,p,f) > \r 1 (E)\ + \f\ + (d-l)\a\-2p + 2, 

with equality if and only if a = /3 € {id —>•/}. 

Counting the number of cycles of /, we have that \f~ (E) \ + \ f\ ^ 2p — 2, with equality iff 
/ = I. We conclude that S(a, f3, f) ^ 0, with equality if and only if a = /3 = id. The result follows 
by plugging the minimizing values into the asymptotic expression (|49[) . d 



Theorem 4.6. In the regime k ~ cn d , wft d € (1,2), f/ie eigenvalues Ai ■•■ A n 2 o/ 
Z = [$ ® $](£„) satisfy: 

• In probability, k\\ — > 1. 

• Almost surely, n i_ 1 ^27=2 ^n 2 Aj converges in distribution to the Dirac mass at 1, Si. 
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5 Asymptotics of the von Neumann entropy 



Using the moment information and the behavior of the lower part of the spectrum deduced in 
Section [4] we now analyze the von Neumann entropy of the random output matrix Z: 

n 2 

H(Z) = -^\ i \og\ i . (51) 

i=l 

In the physical literature, the idea of using conjugate quantum channels to tackle important 
questions, such as the additivity of the minimum output entropy for quantum channels, dates back 
to the work of A. Winter and P. Hay den [3]. To bound the entropy of Z, they use the following 
fact, coming from linear algebra, which is independent of the random model used: the largest 
eigenvalue of Z is larger that the inverse of the dimension of the ancilla space, Ai ^ 1/k. This 
bound (which is actually a bound on the operator norm of Z) yields the following bound on the 
von Neumann entropy: 

H(Z) <21ogA;-^|^ + i. (52) 

One of the main applications of the results in this paper is the fact that our exact spectral 
information yield better upper bounds for the von Neumann entropy in some specific cases. 

Theorem 5.1. The asymptotic von Neumann entropy of the random output matrix 
Z = [$ (g> $](£„) is given by: 

1. If d — (k — c is an integer): 



H(Z) = - 




(53) 

2. Ifde (0,1): 

H(Z) = 21ogfe + o(l). (54) 



Ifd = 1 (see \G$): 



H(Z) 



21ogfe- %■ +o(l) 
21ogn- +o(l) 



if < c < 1, 
if c > 1. 



4. Ifde (1,2): 



H(Z) = 21ogn + o(l). 



(55) 



(56) 



5. Ifd^ 2: 



H(Z) = 21ogn + o(l) 



(57) 



One can make the following very instructive observations about the above theorem. First, 
notice that in all cases where d > 0, the main contribution to the von Neumann entropy is given 
by the lower part of the spectrum, and not by the main eigenvalue. This is in contrast with the 
case of the p-Renyi entropies (p > 0), where the largest eigenvalue gives the main contribution [3]. 

Another important point concerns the fixed dimension case d — and the linear case, d = 1, 
k ~ cn. In these regimes, the entropy defect 21og(fc An) — H(Z) is macroscopic (153|) . (|55[) . This 
improves considerably the naive bound (|52p and may provide more insight into the question of 
additivity of minimal output entropies, as argued in [4] 1211 110] . However, Hastings' techniques 
[4], as well as the recent developments of Aubrun, Szarek and Werner [10] do not seem to apply 
to this linear regime. Studying additivity violations in the linear regime and using the bounds in 
Theorem 15. ll to provide larger violations remain interesting open problems. 
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